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We derive closed-form expressions for 6x6 valence-band and 8x8 Kane Hamiltonians for wurtzite and zinc- 
blende semiconductor heterostructures in cylindrical coordinates, which provide a transparent and computa- 
tionally efficient band formalism for low-dimensional semiconductor systems with the cylindrical symmetry. 
We use a specific combination of the Brillouin zone center Bloch functions as a new basis allowing to treat 
the envelope functions themselves as eigenstates of the operator of projection of total angular momentum on 
the symmetry axis, while the effective Hamiltonians parametrically depend on an eigenvalue of the operator. 
As an application, the developed method is used to map the ground and several excited hole states in an 
isolated wurtzite GaN quantum wire of the hexagonal cross-section to the states in an equivalent quantum 
wire of the circular cross-section. 

PACS numbers: 73.21.-b 



I. INTRODUCTION 

Unique physical properties of artificial electron sys- 
tems with two- and three-dimensional quantum confine- 
ment, realized in semiconductor quantum wires (QWs) 
and dots (QDs), have made them an object of exten- 
sive research, in particular in the last two decades^. The 
current state of the art in the nanostructure technology 
and its steady improvement guarantees that QWs and 
QDs will inevitably become a basis for future genera- 
tions of photonic and electronic devices 2 - - — , with QDs be- 
ing viewed as one of the most promising candidates for 
use in solid-state quantum computation^—. 

Despite the impressive progress in the application of 
atomistic approaches to actual nanostructure o 10-13 , the 
envelope function (EF) method still remains the most 
transparent means to study fundamental properties of 
electron states, and virtually the only efficient and flex- 
ible tool for optimization of nanostructure devices. A 
number of shortcomings of the original metho d i 15 , as- 
sociated with the behavior of the wave functions at the 
heterointerfaces, have been proved to be not fatal when 
recipes of the method's modification have been devel- 
oped and applied 1 ^—. Furthermore, for typical rela- 
tively large systems, with a small ratio of the number of 
heterointerface atoms to the number of atoms located in- 
side homogeneous materials, the interface effects are usu- 
ally weak, and an application of even the classic effective- 



mass approximator 



can be justified 2,0 



25.20 



Already within the EF method, the calculation of elec- 
tron states in QWs and QDs is still an arduous problem 
due to the necessity to solve a set of coupled partial dif- 
ferential equations, supplemented by the strain and elec- 
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tromechanical subtasks22r— . Fortunately, for systems 
that are reasonably well approximated by those possess- 
ing cylindrical or spherical symmetry, one can reduce the 
number of independent variables and greatly facilitate 
the analysis. For a single-band equation, the transforma- 
tion to cylindrical or spherical coordinates is a textbook 
tas k 30 ' 31 . Vahala and Sercel have presented a formalism, 
using the basis of Brillouin zone center Bloch functions 
as eigenstates of total angular momentum, allowing to 
integrate out the angular degree of freedom in the multi- 
band Hamiltonians and leave only the radial, for spher- 
ically symmetric systems, or radial and axial, for cylin- 
drically symmetric systems, variable o 32 ! 33 . Their original 
method, developed for zinc-blende materials in the spher- 
ical band approximatio n 34 ! 35 , has since been used exten- 
sively, with an extension to wurtzite nanostructures 3 ^. 

Following that method, for a cylindrically symmetric 
nanostructure, for example, the effective Hamiltonian for 
each value of projection of total angular momentum on 
the symmetry axis should be obtained separately. As 
a result, the procedure has an infinite set of dissimi- 
lar Hamiltonians as an outcome, with such an intricacy 
reminding the formulation of the Heisenberg's matrix 
mechanics 3 -^. 

We develop an alternative approach to the multiband 
Hamiltonian for radially symmetric systems described 
by 6 x 6 valence-band or 8 x 8 Kane Hamiltonians for 
wurtzite and zinc-blende nanostructures. Considered 
only are systems that have the cylindrical symmetry 
so that the formalism can directly be applied to zinc- 
blende (truncated) conical QDs and cylindrical QWs, us- 
ing the spherical band approximatio n 34 ! 35 , as well as to 
wurtzite cylindrically symmetric nanostructures grown 
along the [0001] direction. We expect that properties 
of corresponding pyramidal QDs can also be qualita- 
tively modeled using the cylindrically symmetric shape 
approximation 3 -!. Our resulting EF equations resem- 



ble the "single-band" Schrodinger equation expressed in 
cylindrical coordinates 3 ^, with the EFs themselves being 
eigenstates of the operator j z of projection of angular 
momentum on the symmetry axis. For each specific case 
(wurtzite or zinc-blende, 6 x 6 or 8 x 8 model), we derive 
a unique effective Hamiltonian parametrically dependent 
on a single parameter that is an eigenvalue of j z . 

In Sec. |ll] we deal with 6x6 Hamiltonians for the 
valence-band states. We start from the Hamiltonian 
in Cartesian coordinates for wurtzite (Subsection III Al) , 
make a cylindrical coordinates transformation (Subsec- 
tion III B|) , include an external magnetic field (Subsec- 
tion QTC]), derive the boundary conditions (BCs) to be 
imposed on the EFs on the symmetry axis (Subsec- 
tion III D|) , and specify what changes in the results are 
to be made considering zinc-blende nanostructures (Sub- 
section III E[) . In Sec. IIIII we generalize the results of 
Sec. |n] for 8x8 Kane model, reviewing the Hamiltonian 
in Cartesian coordinates (Subsection 1111 A[) and present- 
ing the resulting Hamiltonian in cylindrical coordinates 
(Subsection IIII B|) . As a simple application of the ob- 
tained results, in Sec. IIVI we study an isolated wurtzite 
GaN QW grown along the [0001] direction, both with 
the hexagonal cross-section in Cartesian coordinates and 
making the cylindrically symmetric shape approximation 
to use the Hamiltonian in cylindrical coordinates. We 
subsequently compare the spectra of hole states for these 
two cases. Conclusions are given in Sec. |Vj 



II. HAMILTONIAN FOR VALENCE BAND 

A. Hamiltonian for wurtzite in Cartesian coordinates 

We start from the steady state Schrodinger equation 
in the effective-mass approximation 2 ^ 24 ' 38 " — , 



HF 



!HW + H^ + H^ + H( £ M F = EF 



(1) 



constructed using the orthogonal and normalized Bril- 
louin zone center Bloch functions for the potential well 
material, which are written in the matrix form as u = 
(u x u y u z ). They can be chosen real. In Cartesian co- 
ordinates r = (x,y) and z, the functions u x = u x (r, z) 

I 



and u v = u y (r, z) transform as the coordinates x and 
y belonging to the representation Tq of the space group 
Cq v , while the function u z = u z (r,z) transforms as the 
coordinate z (along the c-axis of wurtzite) belonging to 
the representation Ti. The EF matrix F = F (r, z) has 
three components: F = (F x F y F z ) , where T stands for 
the transpose of a matrix or an operator, so that the total 
wave function ^ = ^ (r, z) is as follows: 



*= E u ^ = uF - 



(2) 



The basis functions u are spinless^ 2 -, but each EF element 
Fj is a spinor having two components: 



F 



(«) N 



Fi = J , 



F 



(■■I) 



J = x,y,z. 



(3) 



The Hamiltonian H^ ' entering Eq. (UJ represents the 
position-dependent potential energy of an electron: 



H (°) 



'U v6 

o u v6 o 
u vll 



(4) 



where U v § — U v g (r, z) and U v \ — U v \ (r, z) are the po- 
sition dependent edges of the valence bands T 6 and Ti, 
respectively, which may include an external scalar poten- 
tial. 

The Hamiltonian H^ entering Eq. (fTJ) defines the 
spin-orbit interaction, which is taken into account in the 
first-order perturbation theory, neglecting small terms 
linear in the momentum operator— "—: 
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where a x , a y and a z are the Pauli matrices^: 
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and real A2 = A2 (r, z) and A3 = A3 (r, z) are the pa- 
rameters of the valence-band spin-orbit splitting— "—. 
The Hamiltonian H^ of the Eq. Q is: 



'L 1 k 2 x +M 1 kl + M 2 k 2 z 
HW = [ N x k x k v 

N 2 k x k z + N 3 k x 



Nik x k y N 2 k x k z - N 3 k x 

M x k 2 x + L x k 2 y + M 2 k\ N 2 k y k z - N 3 k y 
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which presents the kinetic energy, where 
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(9) 



with A\ , A%, . . . Aj being real material parameters in con- 
ventional notation s 38 ' 39 ' , m® is the free electron mass, 
and hkj = — iftVj, where j = x,y,z, is the momentum 
operator. Not all parameters are independent: 



Lx-M^Ni, 



(10) 



due to six-fold rotational symmetry of wurtzite 3 - 9 -, which 
actually leads to higher isotropic symmetry of the spec- 
trum in bulk materials 40 . 

A common mistake is to treat the above parameters as 
position-dependent due to variation of the chemical com- 

I 



position of the heterostructure, while being limited with 
the second-order differential equations of the effective- 
mass approximation 19 . Here we assume them as the pa- 
rameters of the potential well material. 

To complete the definition of Eq. (fTJ), the strain Hamil- 
tonian H^ has a structure resembling that of Eq. ([7]) be- 
cause the tensor of deformation potential has the same 
transformation properties as the tensor of the reciprocal 
effective mass entering Eq. (|7|), both being governed by 
the lattice symmetry 3 -. We have: 
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n 2 e xz 



m 2 e z 
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m 2 £z 
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n 2 s xz 
n 2 £ yz 

+ £ vv 



+ he 



2 £zz, 



(11) 



where real parameters Zi, l 2 , mi, m 2 , m, and n 2 are 
expressed via conventiona l 38 ' 39 ' 41 components of the de- 
formation potential tensor as follows: 



h = D 2 + D 4 + D 5 , rni = D 2 + D 4 - D 5 , 
ni=2D 5 , l 2 =Di, m 2 = D 1 + D 3 , 
n 2 = V2D 6 , m,3 = D 2 , 



(12) 



and e = e (r, z) is the strain tensor in Cartesian coordi- 
nates^. Contrary to the effective-mass parameters, 
the deformation potentials can be treated as position- 
dependent because the strain directly modifies the po- 
tential energy of an electron, which is a dominant ef- 
fect in view of the estimation of the terms entering the 
Hamiltonian^ 3 -. For heterostructures composed of related 
materials resulting in a weak position dependence of the 
band parameters, which is a requirement for the effective- 
mass method's applicability, the deformation potentials 
are also expected to have a weak position dependence. 
In this case parameters of the potential well material can 
be used for the whole structure. 

Analogous to the identity of Eq. (ITU)) , an important 
relation takes place: 



h -mi = m, 



(13) 



which secures the cylindrically symmetric strain Hamil- 
tonian for cylindrically symmetric strains. For a cylin- 
drically symmetric system, H", H 1 -^ and e depend on 
the absolute value of the vector r. 



I 

tions u as follows: u = uS, where 

cos(j) sincj) 0\ 

-sin0 cos0 , (14) 

1/ 

which is a standard unitary matrix used to express Carte- 
sian coordinates of a tree-dimensional vector in terms of 
cylindrical ones. Explicit relations of the new functions 
u r and u<p to u x and u y read: 



xu , 



\/x 2 + y 



yu v xu v — vu r 

'' ' u^ = ^JLJL . (15) 



\A 2 + y 2 



Note that u r and u^ are not Bloch functions, in particular 
they are not periodic. 

The total wave function, see Eq. ([2]), is: 



* = uSF = uF, 



(16) 



where we have introduced the modified EF F. The 
effective-mass Eq. (fTJ) now reads: 



HF = EF, 



(17) 



where H = SHS- 1 . 

The potential energy Hamiltonian does not change: 



h( » = sh(°»s- 1 = h(°). 



(18) 



To transform the kinetic energy Hamiltonian presented 
by Eq. (O, we use Eq. (TTUj) and the identities: 



B. Hamiltonian for wurtzite in cylindrical coordinates 

To carry out a transformation from Cartesian (x, y, z) 
to cylindrical coordinates (r,cj>,z), with x = rcos</> and 
y = r sin cj>, it is important to transform all correspond- 
ing elements of the considered system. We thus introduce 
new "basis" functions u = (u r u,f, u z ) related to the func- 



V x = cos(t> V r 



Vy = sin <f> V r 



V,*, 



(19) 



where V,- = d/dr and V0 = d/d<f>. Let us follow some 
algebra necessary to obtain the elements of the Hamilto- 
nian H( fe ) = SH^S" 1 . For example, 



H(i — cos ff}x cos + cos H{ 2 sin + sin H 2l ' cos + sm ^22 sm ^ 
= — L\ (cos Vz + sin V y ) (Va; cos + V„ sin 0) — Mi (cos V„ — sii 



= -Li Vj 



•^ - 



M 1 —Vl-M 2 V 2 



V y sin 0) — Mi (cos V y — sin V x ) (V y cos — V x sin < 

2' 



M 2 V 2 



(20) 
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r«, 



r(fe) 



W, 



rW. 



-ffi 2 = — cos^iJ^ sin0 + cos(/)Hf 2 cos0 — smc^H^i sin0 + sin^-ff^ cos0 

= — L\ (cos V^ + sin V y ) (V v cos — V x sin 0) — Mi (sin V^ — cos V y ) (V x cos + V y sin ( 



ilV r -V,, 

r 



Mi [ -V r V, 



The rest elements are obtained analogously. The full matrix of the kinetic energy Hamiltonian EF fe ) is: 



H« = _ 



^V r 



ii 



^V, 



-MiV r ^ Mi[V? + V r i]+Li- 



m 2 v 2 



\ ^2 [V P + i] V z + iiV 3 [V r + i] 



W 2 ^V, + tJV 3 ^ 



(21) 



M 3 



iV 2 V r V 2 - i7V 3 V r \ 


^V 2 ^V,-z7V 3 ^ 


V 2 + iv r + 4 


+ i 2 V^ 



(22) 



One should bear in mind that V^F = — (V r + r~ 1 ) to 
verify that the Hamiltonian of Eq. (|22[) is Hermitian (note 
also that the parameter N 3 is imaginary). 

In a similar way as above, with the help of the identity 



I 

of Eq. (fT3"]) , the strain Hamiltonian of Eq. dTTI) transforms 
into the Hamiltonian H( £ ) = SH^S -1 : 



f he rr + mie^ + m 2 e zz 
H ■ : ■ [ n\E r <p m\£ r 



n 2 £ r 



n\£rrj> n 2 £ rz 

- l\£^ + rri2£zz n 2 £<p z 

n 2 £^z m 3 (e rr + £4,4,) + hsz 



(23) 



where e — e (r, 0, z) is now the strain tensor in cylindrical 
coordinates^,. We have also used the following identities 
accompanying the coordinate change: 

£xx = £rr cos 2 + e^ sin 2 - £ rtj) sin 20, 
sin + £00 cos 2 + £r<f> sin 20, 

sin 2</> + e r cos 20, (24) 



£ rr 



*x»- 2 

£iez = £rz COS - £0 Z sin 

c- e- oivi ^A _L 

tyz 



e rz sin < 



■ £<t>z 



cos< 



which are in effect for any second-rank tensor. For 
cylindrically symmetric systems, the components of the 
strain tensor, expressed via the displacement vector v = 
(v r ,v^,,v z ), for which v^ = and <9v/<90 = 0, are 42 : 



dv r 

1 / dv r dv 2 
Srz= 2 ( ~dz~ + ~fr 



_ dv z 

£zz — j 

oz 

£r<ft — £(f>z U. 



(25) 



The spin-orbit interaction Hamiltonian of Eq. ([5]) 



I 

transforms into the Hamiltonian U^ = SH^S -1 : 

(0 ~iA 2 a z iA 3 a^ 
iA 2 a z -iA 3 a r ) , (26) 

-iA 3 C0 iA 3 o- r 

where a r = a x exp (i(pa z ) and a,p = a y exp (i<jxj z ). 

The spin-orbit Hamiltonian of Eq. (|26l) depends on 
and does not commute with the operator — iV^. To avoid 
this dependence, we note that 



<J y , 



(27) 



and then construct a unitary transformation of the EF: 
SF = F, (28) 



where 



e l * CT « 
e^° 



(29) 



The effective-mass Eq. (JTTJ) now reads: 
HF = EF, 



(30) 



where H^SHS- 1 . 

The spin-orbit interaction Hamiltonian of Eq. (I26[) re- 
covers its initial form given by Eq. ([5]): 



hW = sh^st 1 = h^ 



(31) 



The potential energy Hamiltonian of Eq. ([18]) and the 
strain Hamiltonian of Eq. f|23[) remain unchanged: 



H(°)+H< £ ) = S (H<°) + H« J S- 1 = H<°) + H«. (32) 
With the identity: 



e^'V^e - **' 



V - -cr 2 



(33) 



the kinetic energy Hamiltonian H (fc ) = SH^S -1 still 
has the form given by Eq. (|22[) where the operator V^ is 
replaced by the matrix operator: 



V,.- 







(34) 



We can express this fact by using the notations: 

H< fc > = H< fc >(V*->V*). (35) 

For cylindrically symmetric systems, the Hamiltonian H 
now commutes with the operator of the z-projection of 
the angular momentum j z = — iV^ (expressed in units of 
K). 



Note, that we could bypass the unitary transformation 
expressed by Eq. (|28|) but use the reciprocal transforma- 
tion of the operator j z : 



J, 



e-*»"*j»e* 



1 



: a z, 



(36) 



explicitly producing the operator of the z-projection of 
total angular momentum J z that commutes with the 
Hamiltonian H. 

Being commuting operators, H and j z have common 

eigenfunctions^i, so that the EF F = F (r, <f>, z) can be 
chosen in the following form: 



^imtfi 



f(r,z), 



(37) 



where m, being an eigenvalue of the operator of the z- 
component of electron's total angular momentum, should 
be half-integer: m = ±1/2, ±3/2, . . . The function f = 
f (r,z); 



f = ( f^ f^ f^ f^ f^ f. 
satisfies the equation: 



(d) 



(38) 



H (°) + H^ + H^ + B$ ) f = Ef, (39) 



with H(°), H( ff ) and H^ being given by Eqs. 0}, © and 
(|23|i. respectively, the latter should be supplemented by 
Eq. (1251) . Having m as a parameter, the kinetic energy 
Hamiltonian is: 



fin [V? + V r i] - iVh^- + M 2 V\ iI*V r f - iM x [f\7 r + $\ 
fiW = i iL x [=V r + ^] - »AfiV r = Mi [V 2 r + V r ±] -L x ^- + M 2 V 2 
N 2 [V r + i] V z + iiV 3 [V r + i] iN 2 fV z -N 3 f 



where the matrix m has a simple form 





_ l 
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m 



M 3 



iV 2 V r V z - iJV 3 V r 
iN 2 fV z +N 3 f 



L 2 V 



(40) 



I 

are not periodic, it is more convenient to use the old 
basis function u, which are periodic and orthonormal, 
U\\ recovering the old EF F: 



The total wave function ^> and the EF F are related to 
each other via the conventional basis functions u or the 
modified ones u as follows: 



(42) 



* = uS^S^F = uS^F, 



where the unitary matrices S and S are given by Eqs. (|14p 
and (|29p . respectively. Since the functions u r and u<j> 



S _1 S _1 F, 



(43) 



as soon as the effective-mass equation has been solved 
and the EF F has been found. 



C. External magnetic field 



where go ~ 2 is the free electron g-factor, \xb is the Bohr 
magnetron, and Q is a material specific parameter—: 



An external magnetic field B z , applied along the z- 
axis, does not break the cylindrical symmetry of the sys- 
tem if such exists. We use the symmetric gauge for the 
vector potential A = (A x ,A yi A z ) = (—yB z /2,xB z /2,0) 
in Cartesian coordinates, or (A r , A^, A z ) = (0,rB z /2,0) 
in cylindrical coordinates. To include the magnetic field, 
we make the following substitution in the Hamiltonian of 
Eq. 0: 



k ] ~> K 3 = k 3 + Y C A ^ 3 = X ' y ' Z ' ( 44 ) 



where — e < is the electron charge, and c is the speed of 
light in vacuum. All non-commutative products should 
be replaced by the symmetrized products 23 : 



K x K y -> -{K x , K y }+ = - (K x K y + K y K x ) . (45) 



ft 2 ^^ r> yx) n vu > 

~ ml*-? E v6 



Pxj Pjy 



E-y 



(47) 



where E v6 is the energy of the valence T 6 state, with 
Ej being the energy of 7 band, for the potential well 

material, p?, are the matrix elements of the momentum 
operator between the Bloch zone-center functions w 7 and 

Uji, where j,j' — x,y, that is pV = (iiji\hkj\u 7 ). 

Let us compare Eq. (J47I) with a similar expression 3 ^ 
for the material parameter A^: 



Ai 



fi 2 J- x )J-y) 1 Jv)J- x ) 

II ^— % Pxj Pjy ~r Pxj Pjy 



m o^ 



EL 



E^ 



(48) 



Jv) 



and note that p x ~' = for the nearest 7 bands, which 
belong to the representation T±. This indicates that the 
approximation Q « N\ may be satisfactory. 

In cylindrical coordinates, for the Hamiltonian of 
Eq. (|2"21 we symmetrize the products of the operators 
Vci and V r : 



In addition, the Hamiltonian H of Eq. (fTJ) must include 
the effective Pauli term H*- 5 ': 



V V r -^ -{V,i,V r }+, 



and make the substitution: 



(49) 




HW = 



v ^ V, 

so that, for example, 



B z r 2 , 
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HW = - V + -— B z r\V r ) -± + _± V + — £ z r z ,V r ^ + -f V + — B z r 



2hc 



Mi 



2hc 



(51) 



The remaining arguments directly follow the ones for 
the case without magnetic field. Thus, the unitary trans- 
formation of Eq. (|28|l brings about the result of Eq. (|35|) . 
The effective Pauli term of Eq. (|46p remains invariant: 



H< B > = SSH^S^S" 1 = H^ 



(52) 



The resulting full Hamiltonian commutes with the oper- 
ator — iV<?i, if the system is cylindrically symmetric. The 
EF has the form of Eq. (|37l) . where the function f is 



obtained by solving the following equation: 

H (°) + H< ff > + gW + H< B > + H%A f = Et, (53) 



where H<°), Bl^ , ET>\ and H( s ) are given by Eqs. g}, 
©, (|23|) . and (|46|) . respectively, while keeping in mind 

Eq. (|25|). The kinetic energy Hamiltonian Hi, is pre- 
sented by the matrix operator: 
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-A 2 V r 



tJV 3 V r 
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I ( ^--Mi[V 2 + V r i]-M 2 V 2 

*A 2 ^V 2 + A 3 ^ M 3 



- A 2 V r V z + iA 3 V r 
^A 2 fV 2 -A 3 ^ 
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(54) 



where the matrix operator m is: 



2hc z ' 







^S z r 2 +m+i 



Also holds Eq. (02 



We will have the same result if we use the continuity 

(u) 

of the function Fy . For the down elements it holds 



(55) analogously: 



/ 1 (d) | r=0 = / 2 (d) | r=0 = 0, \v\?l, 



(63) 



D. Structure of function f (r, z) at r — 

We will determine the BCs imposed on the function 
f (r, z) entering Eq. ([39]) or (|53[) on the symmetry axis 
at r = 0. Having done this, to specify the eigenvalue 
problem, one should also complement Eqs. ([39]) and ([53]) 
with BCs at infinity, e.g., f(r,z) \ r -+oo — > for discrete 
spectrum. 

The BCs at r = can easily be obtained from the 
continuity of the function F, which satisfies Eq. (Q}, and 
its first partial derivatives V^F for j = x,y. We have 
the following explicit expressions connecting the old EF 

F and the modified one F, obtained from Eqs. (|37| and 
dH: 



To obtain the BCs imposed on the derivatives of the 
EF, we use the following continuity conditions for r — > 
and any </>: 

* j V ' ^) \x= r cos (p, y=r sin (j) ' * j \ 5 ^ ) \x=r,y=0-i \P^) 

where j — x,y. For the element /g , using its repre- 
sentation via Fz in Eq. ([56]) as well as the conditions 
of Eq. (|M|) and the identities of Eq. (fTi?]) . we have in 
particular: 



dr r ~ J,i r dr 

as r — > for any </>, from which we conclude that 



9/ W^cos0-^/r^e-^ ^, (65) 



r( M ) Jn<t> , 



v^Fi") = jJ'ieW cos0 - / 2 (u) e 
v/^r>M = /{"W sin0 + / 2 (M) e^ cos./), 



(56) 



27T.F. 



(u) 



f: 



( u ) iv4> 



for the upper elements of the envelope-function spinor 
components, where fi = m— 1/2, and similar relations for 
the down elements, where the EF index (u) is substituted 
by (d) and v = m + 1/2 supersedes fi. 

The continuity of the function F in particular on the 
symmetry axis means that 



F(r,z)| 



x—r cos <p.y—r sin <fi 



->F(r,«)l 



x—r,y—0i 



(57) 



as r — > for any angle cf>. Then, using Eqs. (151)]) and ([5 
for the third component /g we have: 

e^/ 3 (M) (r, *) -> f { 3 u) (r, z) , r -> 0, V 0, (58) 

which is an identity if [i = 0, but definitely 

0, /^0. (59) 



/3 |r=0 



For the down element of the spinor, we obtain analo- 
gously: 



jfi=o = 0, z/^0. 



(60) 



From the continuity condition of Eq. (|57p . using the 
first identity in Eq. (|56p. we have for any 0: 



jfV^ cos - /fV^ sin -» / 1 (u) , r 
which results in the following: 



II \r=0 
p(«) 



f(«)| 

h \r=0 

(«)| 



o, ImI^i, 



/a |r=o = «A*/i |r=o, ImI = 1 



0, (61) 



(62) 



,(«) 



V r ./ 3 w | r =o = 0, ^ = 0. 
Analogously, for the down element: 

V r ./ 3 (d) | r .= = 0, v = 0. 



(66) 



(67) 



For the elements /{ and / 2 , we use the continuity 
conditions of Eq. (|M|) . the first and the second identities 
of Eq. ([56]) along with the identities of Eq. ([19]). We 
obtain in particular that for any </> and r — > 0: 



^e^ cos 2 , 
or 



U J2 J, 

dr 



e ^ sin © cos ( 



Ju) 

^— e 1 ^ (sin 2 (j> - in sin cos < 
r 

/ ; (u) 



(68) 



2 —e 1 ^ (in sin 2 <£ + sin </> cos 0) -> -4f 
r or 



(u) 



For |/x| = 1, using Eq. (1521 . we have: 



(«) 



f(«) 



^e^cos 2 0- £§_e^sin^cos^-). ^, (69) 
or ar or 



which is satisfied for any <\> as r — > if 

V r / 1 (u) | r=o = V r / 2 (u) | r =o = 0, |/i| = l. (70) 
For the down elements of the spinor components: 

V r ./ 1 (d) | r .=o = V T ./ 2 (d) | r=0 = 0, \v\ = \. (71) 



The conditions of Eqs. ([59]) . (j60]l . ([62]) . (j63 |> . (|66]) . (J67]) . 
(f7D|) . and (I7T1) constitute the complete set of the BCs 
that must be imposed on the function f (r, z) entering 
Eqs. ([39]) and ([53]) on the symmetry axis. They are sum- 
marized in Table [J for convenience. 



TABLE I. Boundary conditions imposed on elements of f (r, z) 
at r = for different eigenvalues m of the operator of projec- 
tion of total angular momentum on the symmetry axis. 
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E. Valence-band Hamiltonian for zinc-blende 

The above results can be easily adapted to cylindri- 
cally symmetric nanostructures composed of zinc-blende 
semiconductors. The proper modifications are listed be- 
low. 

The basis functions u Xl u y and u z now belong to the 
representation T 15 of the space group T^, so that in the 
potential energy Hamiltonian of Eq. (Q| we put U v e — 
U v i, while in the spin-orbit interaction Hamiltonians of 
Eqs. © and {26]) we have A 2 = A 3 = A/3, where A is 
the valence-band spin-orbit splitting in cubic materials 22 . 

For the parameters entering the kinetic energy Hamil- 
tonians of Eqs. ©, ([22]). (]40]) , and ([54]) we have the fol- 
lowing identities: 

Lx = L 2 , Afi = M 2 = M 3 , N x = N 2 , N 3 = 0. (72) 

Three parameters L\, M\, and N\ are expressed via the 
Luttinger parameters 71, 72, and 73 as follows^ 3 -: 



— =~(£i + 2M 1 ), ^- = -i(i 1 -M 1 ) 



2too 



73 



1 



(73) 



2?7lo 6 



:Ni 



Equation (|10[) is not an identity for zinc-blende, equiv- 
alent to 72 =73, but the recipe for the spherical band 
approximatio n 34 ' 35 . 

Analogously, in the strain Hamiltonians of Eqs. (fTTj) 
and (12"31 we have: 



h=h 



777,1 = TO2 = TO3, 77 1 = 712- 



(74) 



The valence-band deformation potentials a, b, and d are 
expressed via these parameters as follow s 48 ' 49 : 



1 



1 



(Ii+2mi), b=-(l 



1 



mi) 



V3 



m. (75) 



The validity of Eq. 
alent to 



13]) is not hence guaranteed, equiv- 



d = V3b. 



(76) 



Whether Eq. (f76"|) , which expresses the spherical deforma- 
tion potential approximation, is acceptable or not should 
be analyzed for each specific case. In Table [TTI we present 
a comparison of the deformation potentials ^/5b and d 
along with a comparison of the Luttinger parameters 72 
and 73. It can bee seen that the spherical deformation 
potential approximation may be as good as the spheri- 
cal band approximation, while there exist materials for 
which one or both of these approximations are poor. 



TABLE II. Comparison of the valence-band deformation po- 
tentials y/3b and d, and the Luttinger parameters 72 and 73 
for some zinc-blende semiconductors. The material parame- 
ters are taken from Ref. [50l . 

Parameter GaAs AlAs InAs GaP A1P InP GaN A1N InN 

-2.3 -1.8 -1.6 -1.5 -2.0 -2.2 -1.9 -1.2 

-3.4 -3.6 -4.6 -4.6 -5.0 -3.4 -10 -9.3 

0.82 8.5 0.49 0.71 1.60 0.75 0.47 1.26 

1.42 9.2 2.93 1.23 2.10 1.10 0.85 1.63 

1.17 0.87 0.60 0.56 0.69 1.12 0.33 0.22 

0.58 0.92 0.17 0.58 0.76 0.68 0.55 0.77 



b (eV) 


-2.0 


d(eV) 


-4.8 


72 


2.06 


73 


2.93 


Vsb/d 


0.72 


72/73 


0.70 



III. HAMILTONIAN FOR KANE MODEL 
A. Cartesian coordinates 

The Kane model^ takes into account the direct in- 
teraction of the valence and conduction bands exactly 
while the interaction mediated by remote bands is treated 
to the second order within the Lowdin perturbation 
schem e 38 ' 45 . The model is used for narrow bandgap sys- 
tems when the energies of interest (e.g. band offsets) 
are comparable with the bandgap, but other bands can 
still be treated as remote. There exists a controversy re- 
garding the model for wurtzit o 46 ' 47 , though concerning 
minor contributions yet currently lacking experimental 
values for the corresponding material parameters. Fur- 
thermore, for heterostructures composed of wurtzite type 
semiconductors the potential of variation of the chemical 
composition acquires an "unusual" non-diagonal contri- 
bution due to the fact that the conduction and one of the 
valence band edge Bloch functions belong to the same 
representation Ti of the space group Cq v . We highlight 
all these points below. 

The steady state EF equation of the Kane model can 
be written as follows: 



H'F' = EF' 



(77) 



where 



H' 




F' = 



F c 



(78) 



and H^ is the Hermitian conjugate of H c „. The con- 
duction- and the valence-band EFs are F c = F c (r, z) 
and F v = F v (r, z), respectively, the latter has the com- 
ponents F x , F y , and F z . As in Eq. ([2]), the total wave 
function Vf is: 



* 



j=c,x,y,z 



UjFj 



u'F'. 



(79) 



which explicitly depends also on the conduction-band 
zone-center Bloch function u c = u c (r,z), so that u' = 
(u c u x u y u z ) and F' = (F c F x F y F z ) . As above, the 
basis functions u', which can be chosen real, are spinless 
while each EF element Fj is a spinor: 




J = c,x,y,z. 



(80) 



We now proceed to the description of the matrix 
Hamiltonian H' of Eq. (|78[) . The intra-conduction band 
block is: 



H c =A[k 2 z + A' 2 (k 2 x + kl) + U c 

+ »l£zz + «2 (£xx + £yy) , 



(81) 



J 



where U c = U c (r, z) is the position dependent edge of 
the conduction band Ti, which also includes an exter- 
nal scalar potential if it is present, a\ and a 2 are the 
deformation potentials for the conduction band. The pa- 
rameters A' x and A' 2 are expressed via the components 
1/mi and l/ni2 of the tensor of the reciprocal effective 
mass for the conduction band in the single-band approx- 
imation and the Kane parameters (can be chosen real) 
P\ = —ih(u c \hk z \u z ) /rriQ and P 2 — —ih(u c \hk x \u x )/mo 
as follows: 



A[ 



P? 



2m\ E c — F v \ ' 



A' = 



P 2 

r 2 



2m 2 E c — E v 



(82) 



(83) 



where E c and E v \ are the energies of the conduction Ti 
and valence Ti states, respectively. 
The intra-valence band block is: 



H„ = H<°> + H^ + H< e > + H'<*>, 



(84) 



where H^, U.^ and H^l are given by Eqs. gj), © 
and (fTTj) , respectively. The form of Eq. ([7]) is also pre- 
served for the kinetic energy Hamiltonian H'( ) while 
some of the parameters are adjusted to exclude contri- 
butions from the interaction with the conduction band: 



I L[k 2 x + Khk 2 y + M 2 k 

H'W = J N[k X ky 

N' 2 k x k z + N 3 k x 



1 V -^ ft x fcy 

M x k 2 x + L[k 2 + M 2 k 2 
N' 2 k y k z + N 3 k y 

I 






N 3 k x 

N 3 ky 



(85) 



M 3 (k 2 x 



k l) 



P-l^-z, 



where 

L\ =L X + 

N[=N 1 + 

N'=N 2 + 



P 2 

r 2 


P 2 

T' T - I l 


E c — E v q 
P 2 

r 2 


2 £ ' TP V 
rj c — tiivi 


E c — E v q 
P1P2 {2E C 


— E v \ — E v e) 



For the interband block H 
Eq. (UHl) we have: 



where 



(86) 



2 (E c — E v q) (E c — E v i) 

of the Hamiltonian of 



H OT — I H cv i H, 



cd2 



H, 



cv3 J > 



(87) 



H cvl = iP 2 k x + B 2 k x k z + b 2 e xz , 
- iP 2 k y + B 2 k y k z + b 2 s yz , 



H, 



(89) 



H cv3 = iP x k z + B x k 2 z + B 3 (k 2 x + k 2 y ) 

+ hSzZ + h (e xX + Eyy) + U cz , 



(90) 



where we have neglected strain-induced k • p-interaction 
term o 38 ' , proportional to e, which are very small be- 
cause for typical semiconductor heterostructures e ~ 
10~ 2 . The parameters entering Eqs. ([88|) - (|90| have the 
following meanings: b\ — (u c \D zz \u z ), b 2 — (u c \D xz \u x ) 
and 63 = (u c \D xx \u z ) are the interband matrix elements 
of the deformation potential tensor D 38 i 48 . For the pa- 
rameters B\, B 2 and B 3 we have: 



B 1 



2ml *-> 

u 7 



1 



E v \ — EL E c — E-, 



v (z) v {z) 

fci r^z ' 



(91) 



B, 



H 2 y, 

2m 2 ^ 



E~ 



E r — E~ 



E v \ 



(92) 
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E„-\ — £L 



En — E~ 



P^tfl (93) 



The Kane Hamiltonian given in Ref. |46|, Eq. ( A4) , is 
not applicable to wurtzitc. The terms of the conduction- 
valence interband interaction mediated by remote bands, 
which are proportional to the parameters B\ and B 2 
there, are proper not for wurtzite but zinc-blende 
materials-^, after elimination of a misprint in the ele- 
ments (H 2 )u and (#2)41, and putting B\ =B%, On the 
other hand, the Kane Hamiltonian presented in Ref. |47I . 
Eq. (4), is lacking the term proportional to B\ in our 
Eq. dHOD. 

The interband potential U ez = U cz (r, z) entering 
Eq. ([TO]) can be associated with the difference in the 
Bloch functions for the semiconductors composing the 
heterostructure. We can use U cz — for the space region 
of the reference (potential well) semiconductor, which has 
the periodic lattice potential U\ — U\ (r, z), while for the 
space region of the barrier material with the lattice po- 
tential U% = U 2 (r, z) we hav e 19 ' 24 U cz = (u c \U2 — U\\u z ). 
This step-like contribution exists even for heterostruc- 
tures with graded composition, so that this interband 
offset can be estimated knowing only bulk parameters of 
the semiconductors. It will be done elsewhere. 

If we adopt the cubic approximation^— , which im- 
poses the zinc-blende crystal symmetry on potentials of 
the wurtzite semiconductors, all parameters correspond- 
ing to the controversial terms as well as the interband 
offset will turn to zero: 



matrix 



(l 



S' = 



1 

s. 







o\ 



cos0 sin0 
— sin <fr cos <fi 

\o 01/ 



(95) 



plays the role of the matrix S for the valence-band case, 
see Eq. (fl"4|) . To avoid the (^-dependence in the spin- 
orbit Hamiltonian after the unitary transformation real- 
ized with the matrix S', we also use the unitary transfor- 
mation defined by the unitary matrix S', with the coun- 
terpart of Eq. (|29j) from the valence band case: 



S' 



/Vt CT * \ 

e l i a * 

e l 2 a * 

V e'&'J 



(96) 



The resulting doubly transformed Hamiltonian H' of 
Eq. CZD is: 



H' = S'S'H'S'^S'- 1 = 



H a Hf, 



, ti. cv rl t , 



where 



H c = - A'^ 2 Z - A 2 



9 1 1-2 



+ U C + a x e zz + a 2 {e rr + e^) 



(97) 



(98) 



B\ — B2 — B3 — bi — 62 



U c , 



0. 



(94) 



Also, N 3 — hp z x J /mo — 0, and E v \ = E vG in the de- 
nominators of Eqs. (j8"2")l . (j8"3")l and (I5rj|) . Interesting, that 
contributions of similar origin in zinc-blende materials ■ - , 
which disappear for wurtzite, can hence be treated as 
small on reciprocal grounds, rather than via invoking 
an affinity of the zinc-blende and the diamond lattices. 
In what follows we use the approximations given by 
Eq. ([MD. 

For a zinc-blende nanostructure, we use the relations 
given in Subsection III El along with a\ — a 2 in Eq. (J8T|) , 
mi = ra 2 , P\ = P2, and E v i = E v § in Eqs. 
(1891). and (1901). 



B. Cylindrical coordinates 



The procedure of transformation of the Kane Hamil- 
tonian to cylindrical coordinates closely follows the one 
for the valence band case. We accordingly introduce the 
zone-center functions u' = {u c u r u^ u z ) related to the 
functions u' as follows: u' = u'S', where the unitary 



with 



H C y — I H cvl H cv2 H cvi 



H cvl = P 2 V 



Note that 



H cv2 = Pi-Vs, H cv3 = PiV 2 



H\ vl = ~P 2 V r 



while 



H 



t 

cvl 



-P2-V 6 , II 



cv3 



-PiV 2 



For the intra-valence band block H„ we have: 
I„ = H<"» + H< ff ' + H< 8 > + fi'W, 



(99) 



(100) 



(101) 



(102) 



(103) 



(104) 



where H( ', H^ and H( e ) are given by Eqs. (0), © 
and (|23l) . respectively, bearing in mind Eq. (j2"5")l . For the 
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kinetic energy part H'( fc ) we use the matrix Hamiltonian 
of Eq. (J22J with the substitution V^ -► V^, see Eq. (1351) . 
along with the substitutions L\ — > L' 1; L2 — >• L' 2 N\ — >• 

iV{ and N 2 ^ N 2 , see Eq. 



The external magnetic held B z is treated analogously 
to the valence-band case, via the substitution given by 
Eqs. (|49p and (|50|). and an inclusion of the following Pauli 

term H^ 5 ' in the Hamiltonian H': 



H'W 
















—iLQ B z 

SOME „ D 





2 



(105) 



**£*/ 



which is invariant under^the transformations realized 
with the matrices S' and S'. 

For a cylindrically symmetric system, the Hamiltonian 

H' commutes with the operator — iV^. The EF F' = 
S'S'F', which satisfies the equation: 



H'F' = EF', 
can be chosen as follows: 



F' = 



^im<ft 



?(r,z), 



(106) 



(107) 



where m = ±1/2, ±3/2, . . ., while the eight-element 
function f (r, z): 

f = (7 c (u) rt d) f[ u) i[ d) A u) A d) ft ] 4 d) ) T , (108) 



does not depend on </>. 

The BCs for the valence-band elements of the function 
f on the symmetry axis r = are the same as those 
given in Table U of Subsection III PL The expressions for 
the conduction-band elements of f via the old EF F': 



F (u) 



i{m-l/2)4> i(m+l/2)<t> 

/("), F^ = C = fi d \ (109) 



'27T 



formally coincide with the relations for the valence-band 
elements f^ and /g , see Eq. (|56|) . resulting in the same 
BCs at r = 0, which are presented in Table UTT1 



TABLE III. Boundary conditions for the conduction-band el- 
ements of f ' at r — for different eigenvalues m of the opera- 
tor of projection of total angular momentum on the symmetry 



1/2 



-1/2 3/2 -3/2 other half-integers 



V r / C (u) 

fi d) 



















v r /, 



(d) 



IV. HOLE STATES IN WURTZITE QUANTUM WIRES 

The cylindrical coordinates representation might lead 
to no practical advantage over the Cartesian one if it 
were shown that actual nanostructures, such as wurtzite 
GaN/AIN QDs of pyramidal shape with a hexagonal 
base^i or zinc-blende InAs/GaAs square-base pyramidal 
QDs 52 , could not ever be approximated as cylindrically 
symmetric. Zinc-blende systems are of special concern 
here due to the necessity to concede to the cylindri- 
cal shape, spherical band and spherical deformation po- 
tential approximations, let alone issues 5 ^ not taken into 
account within the classic EF method. On the other 
hand, one may expect that wurtzite-based systems grown 
along the [0001] direction are very close to cylindrically 
symmetric ones due to the cylindrically symmetric band 
structure, in the quadratic effective-mass approximation 
including the Kane model, and the geometric shape rea- 
sons. 

To check this assumption, we present a spectrum of 
the hole states, as a function of the wavenumber k z , in 
an isolated wurtzite GaN QW grown along the [0001] 
direction that has a regular hexagonal cross-section. The 
system itself is of interest^. To calculate the spectrum, 
we have used the 6x6 valence band Hamiltonian H in 
Cartesian coordinates, see Eq. ([IJ. We have also made 
a cylindrically symmetric shape approximation for the 
QW, preserving the same cross-sectional area, see Fig. (TJ 
and have calculated the hole spectrum using Eq. (|3"9"j) in 
cylindrical coordinates. 

We have used the following parameters. The radius 
of the cylindrical QW is R = 2 nm. The same cross- 
sectional area has a hexagonal QW with the hexagon 

side length L = R^2ir/3V3 ~ l.lOOi? = 2.2 nm. The 
electron band structure parameters of wurtzite GaN used 
in the computation are 5 ^: A\ = -7.21, A 2 = -0.44, 
A 3 = 6.68, A 4 = -3.46 A b = -3.40, A 6 = -4.90. For 
simplicity, we have neglected the spin-orbit interaction, 
A2 = A3 = 0, and also the crystal splitting Ai = 
for consistency, so that E v \ — E v q, with the value be- 
ing used as the zero-energy reference point, and A 7 = 0. 
The upper and the down elements of the function f , see 
Eq. (f3"9"]l . are then not coupled, and the states are char- 
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FIG. 1. (Color online) Two studied types of QWs with 
the same cross-sectional area: (a) hexagonal QW with the 
hexagon side length L, (b) cylindrical QW with the radius 

R = L\/3V3/2tt. 
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FIG. 2. (Color online) Spectra of holes in the hexagonal QW 
with the side length L = 2.2 nm, solid line, and the cylindri- 
cal QW with the radius R = 2 nm, dot symbols, the latter 
characterized by the projection l z of orbital momentum on 
the z-axis. 



acterized by the projection of orbital momentum on the 
z-axis l Zl with l z — m— 1/2 for the upper spinor elements 
of f , and l z = m + 1/2 for the down ones. 

Due to properties of the space group Cq v , all states in 
the hexagonal QW are either simple or have the two-fold 
character—, not considering the spin degeneracy. The 
simple states are expected to correspond to the ones that 
have l z = within the cylindrical QW model. The two- 
fold degenerate states in the hexagonal QW are then to 
be related to the states l z = ±n, for n = 1,2,..., in the 
cylindrical QW. Indeed, as can be seen in Table IIVI and 
Fig. [3J the ground and several excited states in QWs of 
these two geometries not only qualitatively correspond to 
each other, but are also fairly close in energy. An excep- 
tion holds for the states with \l z \ = 3, which correspond 
to the simple states in the hexagonal QW. 



TABLE IV. Spectrum E„ of several hole states (each dou- 
bly spin-degenerate) in a wurtzite GaN QW with the regular 
hexagonal (side length: 2.2 nm) cross-section and correspond- 
ing energies E^ c ' of hole states in the QW with the circular 
(radius: 2 nm) cross-section for two values of the wavenumber 
k z , presented in Fig. [5] The spectrum E^ c ' is characterized 
by the projection Z 2 of angular momentum on the symmetry 
axis indicated. 
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a h = 0. 

b ll.l = l. 

c |J,|=2. 
d |k| = 3. 
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Degenerate with the E^„' level 



V. CONCLUSIONS 

Having composed "cylindrical" combinations of the 
Brillouin zone center Bloch functions to be used as new 
basis functions, we have made a cylindrical coordinate 
transformation of 6 x 6 valence-band and 8x8 Kane 
Hamiltonians for wurtzite and zinc-blende semiconduc- 
tor heterostructures. As a result, for cylindrically sym- 
metric structures (using the spherical band approxima- 
tion for zinc-blende), we have guaranteed that the EFs 
themselves are eigenstates of the operator of projection 
of total angular momentum on the symmetry axis, while 
the effective Hamiltonians parametrically depend on an 
eigenvalue m = ±1/2, ±3/2, ... of the operator. In the 
Hamiltonians, we have taken into account the deforma- 
tion effects and have made an allowance for an external 
magnetic field applied along the symmetry axis. The ad- 
vantage of the presented method over the one developed 
by Vahala and Serce l 32 ' 33 is the formulation of a univer- 
sal Hamiltonian (for each specific case, wurtzite or zinc- 
blende, 6 x 6 or 8 x 8 model), valid for any half-integer 

TO. 

The results have been applied to the problem of hole 
states in an isolated wurtzite GaN QW grown along the 
[0001] direction with a regular hexagonal cross-section. 
Making the cylindrically symmetric shape approximation 
for this QW, we have shown a fairly good correspondence 
of the spectrum in the cylindrical QW to that in the 
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hexagonal one. 
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